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Abstract. Polarizable randomly charged dielectric objects have been recently shown to exhibit long-range 
lateral and normal interaction forces even when they are effectively net neutral. These forces stem from an 
interplay between the quenched statistics of random charges and the induced dielectric image charges. This 
type of interaction has recently been evoked to interpret measurements of Casimir forces m vacuo, where 
a precise analysis of such disorder-induced effects appears to be necessary. Here we consider the torque 
acting on a randomly charged dielectric surface (or a sphere) mounted on a central axle next to another 
randomly charged surface and show that although the resultant mean torque is zero, its sample-to-sample 
fluctuation exhibits a long-range behavior with the separation distance between the juxtaposed surfaces 
and that, in particular, its root-mean-square value scales with the total area of the surfaces. Therefore, the 
disorder-induced torque between two randomly charged surfaces is expected to be much more pronounced 
than the disorder-induced lateral force and may provide an effective way to determine possible disorder 
effects in experiments, in a manner that is independent of the usual normal force measurement. 

PACS. 05. 40. -a Fluctuation phenomena, random processes, noise, and Brownian motion - 34.20.Gj In- 
termolecular and atom-molecule potentials and forces - 03.50.De Classical electromagnetism 



1 Introduction 

The surfaces of dielectrics, crystalline solids and metals of- 
ten exhibit random monopolar charge distributions [Tl[31 
[2] . Random charges may result from adsorption of contam- 
inants and/or the presence of impurities that can generate 
surface charges which depend strongly on the method of 
preparation of the samples. Surface charge disorder may 
also originate from the variation of the local crystallo- 
graphic axes of the exposed surface of a clean polycrys- 
talline sample which induces a variation of the local sur- 
face potential [2111112] ■ The heterogeneous structure of the 
charge disorder on dielectric surfaces and thus its statis- 
tical properties can be determined directly from Kelvin 
force microscopy measurements •[1 . Randomly charged sur- 
faces are equally abundant in colloidal and soft matter sys- 
tems [iliniEj, examples arise in surfactant coated surfaces 
[H] , unstructured proteins and random polyelectrolytes and 
polyampholytes [9]. 

e-mail: a . na j i@damtp . cam . ac. uk (corresponding author) 



A number of authors [2] have pointed out that dis- 
order effects may significantly influence the measurement 
of the Casimir-van der Waals (vdW) forces between solid 
surfaces in vacuum. These forces act between all objects 
and are relatively short-ranged in nature. Recent ultrahigh 
sensitivity experiments of the Casimir force have how- 
ever revealed a residual long-range interaction force which 
dominates at sufficiently large separations, and it has been 
suggested that it is due to surface disorder effects [2|. 

Recently it has been proposed that quenched random 
charge disorder on surfaces as well as in the bulk of dielec- 
tric slabs can lead to long-range interactions even when 
the surfaces are net-neutral [TUl[TT| . These long-range in- 
teractions stem from a subtle interplay between the quenched 
statistics of surface or bulk charges and the image charge 
effects generated by the dielectric discontinuities present 
at the bounding surfaces. 

It was subsequently demonstrated |12) that two ran- 
domly charged surfaces can interact with both random 
normal forces, whose mean value turns out to be finite 
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and long-ranged as noted above, and random lateral forces, 
which-for two juxtaposed planar surfaces carrying statis- 
tically homogeneous random charges-turn out to be zero 
on the average. Both quantities however show sample-to- 
sample fluctuations whose root-mean-square value also ex- 
hibits a long-range behavior with the separation distance 
between the surfaces. 

In this study, we pursue our analysis of disorder effects 
on long-range interactions and investigate the torque act- 
ing on a randomly charged planar or (spherical) dielec- 
tric object mounted on a central axle next to a randomly 
charged dielectric substrate. We show that although the 
resultant mean torque in this system is zero, its sample- 
to-sample fluctuation exhibits a long-range behavior with 
the separation distance. Even more importantly, the root- 
mean-square value of the torque fluctuations scales with 
the total area of the surfaces and thus represents an ex- 
tensive quantity. Therefore, the disorder-induced torque 
between two randomly charged surfaces is expected to be 
much more pronounced than the disorder-induced lateral 
force and may present a more effective method to quantify 
charge disorder effects in experiments. Our results on dis- 
order induced torque may also be related to the so-called 
lock and key phenomena, which underpin highly specific 
interactions between complex biological molecules such as 
proteins, where long-range electrostatic interactions can 
induce pre-alignment which enables complex molecules to 
interact in a biologically useful manner }13j . 



2 Two plane-parallel dielectric slabs 

Consider two plane-parallel dielectric half-spaces with the 
bounding surfaces separated in the z direction by a dis- 
tance I. The surface at z = belongs to a dielectric half- 
space with the dielectric constant £2 and the surface at 
z = I belongs to a dielectric half-space with the dielec- 
tric constant ei (see Fig. [T]). We call these surfaces 5*2 
and Si, respectively. We denote by the dielectric con- 
stant of the intervening material. Let each surface (la- 
belled by a, /3 = 1, 2) have a random surface charge den- 
sity Pa(x) — Pai'^, z) with zero mean (i.e., the surfaces are 
net-neutral) and the correlation function in the plane of 
the slabs (r,r' e Si,S2) 

{Pa{r, z)pfi{v' ,z')) = 6ap gasS{z ~ la)S{z' - If}) Ca{r - r'), 

(1) 

where we define I2 — and = L In addition, we assume 
that the charge distribution on surface 5*1 is restricted to a 
finite area A. The surface Si is assumed to be mounted on 
an axle that allows for rotation around its central symme- 
try axis. In the case where the random charge is made up 
of point charges of signs ±e of surface density rias , we may 
write the variance of the charge disorder 
The correlation function C(r — r') has dimensions of in- 
verse length squared, meaning that its two dimensional 
Fourier transform is dimensionless. Typically the values 
of Tis for reasonably pure samples are smaller than the 
bulk disorder variance which is usually in the range of 




Fig. 1. (Color online) Schematic representation of two dielec- 
tric half-spaces at a separation distance of / carrying quenched 
random charges (shown by small dark and light brown patches) 
on their bounding surfaces. The charges on the lower surface 
are assumed to be distributed in a circular region of area A, 
which is then taken to infinity. The lower surface is allowed to 
rotate around the central z axis as shown in the figure. 

between 10"" to 10"'^ nm^ (corresponding to impurity 
charge densities of 10^" to lO^^ e/cm^ pUfTSlfTOlITT] ) . 
The electrostatic energy of the system is given by 

E=^Jdxcj,{x)p{^) (2) 

where /9(x) is the total charge density and 

0(x)= J dyG(x,y)/,(y) (3) 

is the electrostatic potential, while G(x, y) is the Green's 
function obeying 

£oV-[£(x)VG(x,y)] = -<5(x-y) (4) 

with £(x) being the local dielectric function. Upon chang- 
ing the charge distribution the corresponding change in 
the energy of the system is thus given by 

SE^ I d^dySp{^)G{^,y)p{y). (5) 

If the charge distribution on the surface Si , pi (x) , is made 
up of point charges, we have 

PiW = ^ 9„(5(x-x„), (6) 

neSi 

where g„ is the charge at the site x„. Now on rotating 
the surface Si by an angle 6 around its symmetry axis, 
that is to say in the direction perpendicular to the normal 
between the bounding surfaces of the two dielectric media, 
we find that the new charge distribution is given by 

P'l W = ^ ^n)S{z - Z„), (7) 

neSi 

where Rg is the two-dimensional rotation matrix. For an 
infinitesimal rotation angle 69, one has Rse = 1 — 1{69) a2. 
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where 02 



-I 

1 



is the Pauh matrix. This means that 



we can write (assuming the summation over the in-plane 
Cartesian components i,j = 1,2 and using the fact that 
the diagonal elements of cf2 are zero) 



In order to proceed we must assume that the correlation 
length of the charge disorder is much smaller than the 
linear dimensions of the area A on Si which is covered 
by random charges. We introduce the relative coordinate 
u' = r' — s' to obtain 



(5p(x) = (5/91 (x) = {a2)ij 
d 



Jjpim{^) = - dr'du' [rjir'p-u'j- 



d 



Sixs; 



d d 



-Ci(u')]e- 



Zn) = 1(69) {a2)ijrj—pi{Ti). (8) 
or, 



(16) 

where S[ is the shifted integration region over u'. Now 
with the assumption that the correlation length of the 
charge disorder is much smaller than the linear dimen- 
sions of the area A on Si we see that the first term above 



As the surface Si is rotated the self-interaction be- 
tween the charges on each surface is unchanged, thus the 
energy change is only given by the interaction of the charges dominates for a large system. Furthermore, we can take 
and image charges in 5*1 with those in 6*2. We may thus 
write 



the integration over u' to be over to obtain 



SE = i{de) (a2)„- J dr'drdzdz' [r'^-^Piir' , z')] 

xG(r~r';z,z')P2(r,z), (9) 

where r' and r are again the two-dimensional coordinates 
in the planes of 5*1 and S2 respectively and z and z' are 
the respective coordinates normal to the planes. We thus 
note that the integration over the coordinate r' is over a fi- 
nite area A, while that over r is unrestricted. The torque r 
acting on the surface Si is thus given by SE = —{66) t. As 
the charge distribution on the surfaces Si and S2 are un- 
correlated we find that {SE) = —{S6){t) — 0, where (• • • ) 
denotes the ensemble average over the random charge dis- 
tributions. Thus the mean torque is zero. The variance of 
the torque is obtained using ({SE)'^) — {S6)^{t^) as 

see equation (1(J^ 

where as before r', s' £ Si, r,s £ S2 and the summation is 
again over the in-plane Cartesian components i,j,m,p — 
1,2. Again as the charge distributions on the two surfaces 
are assumed to be uncorrelated, the only nonzero correla- 
tions in the above are given by 



kikjyiCi {\^)Kjp 



where 



Kjp = I dr r,r„ 



(17) 



(18) 



This then gives the general result 
dk 



= -9is92 



G(k;0,OG(-k;0,0 



(2^)2 

xC2(k)(7i(k)[(a2) ij{^2)mpkikmKjp\ ■ 



(19) 



This formula can then be applied to a general finite two- 
dimensional charged area A of arbitrary shape assuming 
that it is sufficiently large. An alternative, more geometric 
form of this result is 



dk 

J2ny 



G(k;0,Z)G(-k;0,0 



xG2(k)Gi(k) / dr|kx r|2,(20) 



(P2(r,z)p2(s,C)> 
(pi(r',z')Pi(s',C')> 

This then yields 



g2sS{z)S{0C2{r~s), (11) 
<?i,5(z'-0<5(C'-0Ci(r'-s').(12) 



9isg2s 



Si 

where the integration over Si is about the axis of rotation. 
Specializing to the case where the area A on the surface 
Si is a disc and the axis of rotation is at its center, we 
find that 

ttR'^ . A^ 



see equation (13) 



K ~ S- 

11.Jp Ujp 



= S. 



'47r' 



(21) 



We now write the above result in terms of the two-dimensionaF^^^e ^ 'I disc radius and A = nR is its area. Now 



in-plane Fourier transforms G and G2 of the functions G 
and G2 and carry out the integrations over the coordinates 
r and s of the infinite plane S2 to obtain 



using the fact that {d'2)ij{d'2)mj ~ ~Sim and that G(k) 



(t^) = -(0-2)ij(0'2)mp5ls52s 



dk 



and GQ(k) are functions of |k 

, 2\ A'^9ls92. 



2G(k;0,0 



87r2 



k only, we may write 
dkk^[G{k;0,l)]^C2{k)Ci{k), (22) 



xG(-k; 0,0^2 (-k) J,p„„(k). (14) 

Now we have to evaluate the integral over coordinates of 
5*1 which has the form 



which is the final expression of the torque fluctuations 
in the plane-parallel geometry. This result may then be 
compared with the lateral force fluctuations which were 
derived in the form [T^ 



J jpiin (k) ^ 



5ixSi ^ ^dr. 



d d 



ds 



Ci{v'~s')]e 



2k- (r' — s') 



(15) 



{F. 



ASijgis92! 
27r 



kdk hkj[G{k;0,l)]^Ci{k)C2{k), 

(23) 
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(r^) = ~{a2)ij(a2)mp(^ J dr' dr dz dz' ds' dsd( d(' [r'j -^pi{r' , z')]G{r - r' ; z, z')p2{r, z) 

x[47j?-Pi(s',C')]G(s-s';C,C')p2(s,Oy (10) 

/d d 
dr' dvds dsG{r - r'; 0, l)G{s - s'; 0, l)C2{r - s) Ir'^s' ——Ci{r' - s')l . (13) 
dr'- ds'jn 



leading to an intuitively clear physical relation between 
the lateral force fluctuations and the torque fluctuations 
of the form 



l/^(i)j7(i)u. 



(24) 



where we have used the summation convention and the 
moment of inertia tensor definition as 



dr {S, 



(25) 



The torque fluctuations are thus connected with the lat- 
eral force fluctuations through a geometric factor encoded 
by moment of inertia tensor. This result is completely gen- 
eral and valid for the assumed plane-parallel arrangement 
of the two disorder-carrying dielectric surfaces. 




(26) 



2.1 General scaling of torque fluctuations 

We can now proceed to a general analysis of the torque 
fluctuations in terms of the area of the interacting surfaces 
A as well as the normal separation between them /. As the 
lateral force fluctuations have the following scaling with 
respect to A and I [12] 

A 

it follows that typical torque fluctuations (or its room- 
mean-square Tims) scale as 

Trms-y, (27) 

which is thus extensive in terms of the area A. In other 
words, the magnitude of the variance of the random torque 
exhibits a non-extensive scaling with area. This is as one 
might expect, because torque is determined by the geom- 
etry of the area A even in the limit of large area (which 
is not the case for the random lateral force |12j). How- 
ever the geometry dependence of the torque fluctuations 
and the scaling with the area A are obtained simply from 
the moment of inertia tensor. For instance, if instead of 
a disc-shaped area, one considers a square of area A on 
the surface 5*1 which can rotate about its center, then one 
finds that Kjp = SjpA'^ / 12 and consequently 



24tt 



dkk^[G{k;0,l)]^C2{k)Ci{k). (28) 



The torque fluctuations for a square are thus slightly larger 
than that for a disc of the same area as one would intu- 
itively expect. 



Fig. 2. (Color online) Schematic picture of a sphere {S2) ap- 
posed to a planar substrate (Si) at a minimum separation of I. 
Both objects carry random charges on their bounding surfaces 
(shown by small dark and light brown patches). The sphere is 
allowed to rotate around the central z axis. 

2.2 Torque fluctuations in a homogeneous system 

In the case where there are no dielectric jumps at the 
boundaries (ei = £2 = £m) and the charge disorder is 
uncorrelated Ca{r — r') = (5(r — r'), the expression (22) 
reduces to 

, 2\ gisg2sA'^ 



(29) 



This result can be obtained in a direct, more geomet- 
ric manner, by starting from the variation in the bare 
Coulomb energy upon rotating the surface Si by a small 
angle around its central axis. 



5E 



1 



dx dy 



6pi{x.)p2{y) 



and showing after some manipulations that 



(30) 



(47r£o£„ 



ds dr 



[(r-s)2 + P]3 



(31) 



3 Torque fluctuations in the sphere-plane 
geometry 

Let us now consider a sphere (iS'2) of radius R apposed 
to a planar substrate (Si) at a minimum separation of 
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I (see Fig.[2|, both carrying quenched random charges of 
density (variance) gis and (72s over their surfaces only. The 
substrate is assumed to be fixed and the sphere is assumed 
to be mounted on an axle that ahows for rotation around 
the 2; axis. 

We shall first focus on the case where there are no di- 
electric inhomogeneities at the boundaries (ei — £2 = £m) 
and that the charge disorder is uncorrelated. The more 
general case of inhomogeneous dielectrics will be consid- 
ered later (see below). In the absence of image charges, 
the sample-to-sample torque fluctuations in this system 
can be evaluated exactly as follows. 

The variation in the charge distribution of the sphere 
due a small rotation Sip around its central z axis can be 
written as 5p2{'x) ~ —Sip-^p2(pi), where x — (x, 9, ip) and 



P2 



neS2 



9n 1 2 

sm ( 



S{x - Xn) S{0 - 6n)S{ip - ipn) 



is the sphere charge distribution made up of point charges 
located at positions x„ . The variance of the energy of the 
system thus follows as 



{{SE f) = J dxdx'dydy 



,(pi(x)pi(x'))(<5p2(y)<5p2(y')> 



where for uncorrelated charge disorder assumed here 



(33) 



(pi(x)pi(x')) = gisS{r - r')S{z - I ~ R)S{z' - I - R) (34) 
{Sp2{y)5p2iy')) = {5^f (]^^) ^iv - my' R) 

-6{V~V'). (35) 



.i?2 sine 



dipdtp' 



Here we have introduced x = (r, z) G 5*1 and y = {y, 6, ip) £ 
82- Then using {{SE)"^) = {S(p)'^{T^), we find an expression 
resembling Eq. ( [sl] ), i.e. 



9ls92s 



(47r£o£m)^ 



dr R" sin 6 d9 dip 



[{r - sy + {I + R{1 ~ cosO))^] 



(36) 



where s = (i?sin^ cos 93, i?sin sin (^), or explicitly, via a 
parameterization of r as r = (rcosa, rsina), we find 



see equation (31) 
This equation can be written as 



2TrR^gisg2s 



f 



where the dimcnsionless function f(v) reads 
see equation Tsl 



(38) 



Note that the same expression can be obtained if one 
assumes that the sphere S2 is fixed and the substrate 5*1 is 
allowed to freely rotate around the z axis passing through 
the center of the sphere. 



10' 



10' 



^ 10 
10 



10' 



10' 10' 10 ' 10^ 10 ' 10' 10' 

V 

Fig. 3. (Color online) The log-log plot of the dimension- 
less function /(«) as a function of the dimensionless sphere- 
substrate separation v — l/R (solid line). The dashed (green) 



(32) liii6 shows the limiting behavior f{v) 



for large v. 



In Fig. [3] we plot the dimensionless function f{v) as a 
function oi v — l/R. For large separation or small sphere 
radius u = Z/i? S> 1, the torque variance tends to zero and 
we find that f{v) ~ (1.0470 . . .)/v'^; hence 



9ls92sR 

(47reoe™)2Z2 



(40) 



whereas for small separation or large sphere radius v — 
l/R < 1, we find f{v) - const. ~ 27r(16.1705 . . .); hence 



gisg2sR^ 



(41) 



The dielectrically inhomogeneous case where the di- 
electric constants of the sphere and the substrate are in 
general different from Em and from each other is not tractable 
analytically. It is nevertheless possible to provide an esti- 
mate for the torque fluctuations variance in analogy with 
the results presented in Ref. [12] for the lateral force fluc- 
tuations in the dielectric sphere-substrate system. There it 
was shown, using the proximity force arguments, that the 
dielectric effects lead only to a correction of the prefactor 
of the result obtained for a homogeneous sphere-substrate 
system, where the dielectric-dependent prefactor is exactly 
the same as obtained for two semi-infinite planar slabs 

m- 

This result then suggests the following heuristic esti- 
mate for the torque variance in the dielectrically inhomo- 
geneous sphere-substrate system 



'^R^9ls92s 

-nei 



4|ln(l-Z\iZ\2 



(em +ei)2(£™ + £2)2^1/^2 

(42) 

for OL — 1,2. The above result obvi- 



where — |^ 
ously reduces to the one in Eq. (|37|) when ei = £2 = £ 
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(47reo£r- ''^ 



oc riv r27T 

dr de dip 



m.) JO Jo JO 



R'^r^ sin^eisinV 



_^ ^2 _^ _^ _ 2ri?sin0cos V? - 2R(l + R) cosOY 



(37) 



/(.) . Hdu rde rd^ n^sinW^^ 

Jo Jo Jo [1 + u"^ + {v + 1)2 - 2u sine COS ip-2{v + l) COS ef 



(39) 



4 Discussion and Conclusions 

Generalizing our previous results on the force fluctua- 
tions between charge-disordered interfaces, we have now 
calculated fluctuations in the torque between two coax- 
ial randomly charged surfaces. Surprisingly, the disorder- 
induced torque fluctuations scale differently with the sur- 
face area of the interfaces bearing charge disorder, with 
the root-mean-square torque being extensive in the sur- 
face area A. In our opinion this opens up a feasible way 
to measure charge disorder-induced interactions between 
randomly charged media in a way which is independent of 
normal force measurements and with a higher signal than 
lateral force measurements. 

The measurements we have in mind would of course 
have to be sample-to-sample torque fluctuations. While 
they decay with the separation between the bounding sur- 
faces carrying the disordered charges in a similar way as 
the forces, the magnitude of torque fluctuations is VA 
times larger than the corresponding lateral or normal force 
fluctuations which are always comparable in magnitude. 
Results obtained for the sphere-plane geometry in fact in- 
dicate that the disorder effect on the torque fluctuations 
should be highest in this case, asymptotically leveling off 
at a finite value for small separations. 

That vdW interactions between anisotropic media lead 
in general to torques induced by electromagnetic field fluc- 
tuations was first realized by Weiss and Parsegian [Tni 
[17]. The effect persists not only in the non- retarded but 
also in the retarded limit and should be experimen- 
tally detectable with modern instrumentation [Hj such as 
the torsion-balance-based setups ^20] . It is worth stressing 
that in the case of vdW torques the fluctuations are in the 
field mediating the interactions, while the boundaries do 
not fiuctuate and are not disordered. In the case analyzed 
here however, the effect is very different: the field does not 
fluctuate but its sources indeed are quenched in a statis- 
tically disordered state. The other important difference is 
that vdW torques arise due to the non-isotropic dielec- 
tric response of the interacting bodies, while the disorder- 
generated torque fluctuations are present even between 
bodies with a completely isotropic-but disordered-charge 
distribution, being in this sense obviously more universal. 

In order to detect sample-to-sample variation in the 
torque or the force itself, one would have to perform many 
experiments measuring the value of the torque as one ob- 
ject is gradually rotated above the other. If the measure- 
ment of the torque at any given angle can be done to an 
accuracy greater than the value of the torque fluctuations 
estimated here, then valuable information about the effect 
of charge disorder induced interactions could be extracted. 
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